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TWO-CENTER PROBLEM ORBITS 
AS INTERMEDIATE ORBITS FOR 
THE RESTRICTED THREE-BODY PROBLEM 

by 

R. G. Langebartel 
Goddard Space Flight Center 

SUMMARY 

The restricted three-body problem is treated by keep- 
ing for the intermediate orbit Hamiltonian all terms except 
those representing a generalization of the centrifugal and 
Coriolis forces. The intermediate orbit, therefore, takes 
into account the gravitational forces of both finite masses. 
An approximation is presented for representing a class of 
intermediate orbits similar to the orbits suggested for the 
Apollo project. The perturbation equations are given in 
outline. 
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TWO-CENTER PROBLEM ORBITS 
AS INTERMEDIATE ORBITS FOR 
THE RESTRICTED THREE-BODY PROBLEM 

by 

R. G. Langebartel 
Goddard Space Flight Center 


INTRODUCTION 

The early masters in the field of celestial mechanics, particularly Euler, discovered that the 
problem of the motion of a particle about two fixed mass points ("two fixed centra tT )can be reduced 
to quadratures, although these integrals do not represent elementary functions (Reference 1, p. 
145). Over a hundred years later Charlier (Reference 2, Vol. n) suggested that the two-center 
orbits could be used as intermediate orbits in the restricted three-body problem, and he sketched 
a broad outline for such a procedure- His program was carried out in some detail by Samter (Ref- 
erence 3). He found the calculations lengthy and the perturbations relatively large, which seemed 
to dampen his enthusiasm somewhat for Charlier f s idea. Samter treated the case of the two given 
masses’ being nearly equal, basing his approximations on this circumstance. In the present article 
we investigate the efficacy of Charlier’s method for orbits of the type suggested for the Apollo 
mission. The near- ejection type orbits admit of an approximate representation for the intermedi- 
ate orbits in terms of elementary functions. 


THE EQUATIONS OF MOTION 


Let the two finite masses in the restricted three-body problem be indicated by // and 1 - u. . 
These two masses are assumed to be describing circular orbits about their common center of 
gravity. Let (x, y) bebarycentric coordinates for a system rotating with the masses so that the co- 
ordinates of the particle with mass ^ are (1 - ij., 0) and of the particle with mass 1 - m are (-/ j. , 0). 
The Lagrangian and force function for the restricted three-body problem in this rotating coordi- 
nate system then are (Reference 1, p. 350): 


L = \ (x 2 + y 2 ) + (xy-yx) + U(x, y) , 


(x 2 +y 2 ) 


1 ~ M 


|(x +/x) 2 +y 2 | 1/2 |(x - 1 +m) 2 +y 2 l 1/2 


1 



Since the two-center problem is separable in either elliptic or bipolar coordinates and since 
we shall be using (certain generalizations of) two-center orbits as intermediate orbits we deter- 
mine the Lagrangian, force function, and Hamiltonian in these new coordinates. 

Bipolar coordinates, r x , r 2 , are the distances of (x, y) from the points (-/x, 0) and (1 - 0), 

respectively. Thus 


r ! = |(x + y 2 | 1/2 , 

r 2 = l(x- 1 + /x ) 2 + y 2 | 1/2 

Elliptic coordinates, £ , 77, have simple representations in terms of r x , r 2 : 

cos £ a r x - r 2 
cosh T) - r x + r 2 

Consequently, the transformation formulae from the barycentric coordinates to the elliptic are 




x - -j 


fj. + ^ + \ cos £ cosh 77 


2 sin £ sinh 77 . 


The restricted three-body problem Hamiltonian in bipolar coordinates has the form 


l f , 1 - ^ ~ r/ , ^r j y[(r, + r>-lJ [l - (r, - r^] 

H ' 2 K " 2r x r 2 K ^2 + 7 ^2 + 2r x r 2 ' W 


[1 -M) r i y[(r 1 +r 2 ) 2 -l] [l - (r,-r 2 ) 2 ] ^ ^ _ 


^ r l r 2 


2 r . 


i£ _J£_ , (1) 


where are the momenta conjugate to r x . And in the case of elliptic coordinates with conju- 
gate to £ and fc 2 conjugate to 77 the restricted three-body Hamiltonian is 




t 2 


cosh z 77 


( 2 ) 


2 
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The formulae for the two-center problem are 

L< 2 ) = ^ (x 2 + y 2 ) + UC 2 > = g (cosh 2 77 - cos 2 g )(^ 2 + V 2 ) + U ^ 2 ^ , 

U< 2 > = 


1 “ M 


1/2 


H< 2 > 


|(x +^) 2 + y 2 

2 

cosh 2 r) - cos 2 g 1 1 cosh 2 77 - cos 2 £ 


K 2 + 



2(1 -M) 

+ 2 i X - 


, I 1/2 
+ y 2 

cosh 77 + cos g 

cosh 77 - 

cos £ 

i. 2 

2(1 -m) 



2 <r l a '2 

cos^ 

cosh 77 + cos £ 

cosh 77 - 

cos £ 


(3a) 

(3b) 

(3c) 


Actually, of course, Equations 3a, 3b, and 3c are not strictly the formulae for the two fixed centers 
problem if we adhere to our definition of x, y since these represent coordinates in a rotating sys- 
tem whereas in the Newtonian two-center problem the two masses are fixed in an inertial frame. 
A more nearly precise statement would be that h c2) represents the H of (2) with the Coriolis and 
centrifugal force terms (the terms linear in fc £ ) suppressed. Note, though, that the represent 
different quantities inH andH< 2) since the Lagrangian for the restricted three-body problem has 
terms in £ and v as well as £ 2 and ^ 2 . 


PARTITION AND SEPARATION OF VARIABLES 

In view of (2) we see that Hamilton’s Equation for the restricted three-body problem in ellip- 
tic coordinates is 

aw . 9 \ , f /aw\ 2 /aw\ 2 r n aw 

at + (cosh 2 77- cos <f) + 2 {d v ) - sinh 77 [cosh 77 + (1 - 2/z) cos^J 

- sin £ [( 1 - 2 fi) cosh 77 + cos ^ - 2 jcosh 77 - ( 1 - 2 /jl) cos = 0 . (4) 

The variables <f, 77 do not separate in this partial differential equation although they do in Hamilton’s 

Equation formed from h (2) . Thus, a separable problem for an intermediate orbit would result if 
H be partitioned into H 0 with H 0 = h (2) . However, a larger class of intermediate orbits is ob- 
tained from the following setup: 

"H - H 0 + Hj , 

H 0 = 2 (cosh 2 77 - cos 2 <f) -1 (t>j 2 + fc 2 2 ) - (cosh 2 77 - cos 2 <f)" 1 [u(£) V 1 + v{r/) fc 2 ] 

2(1 -M)_ 
cosh 77 + cos g 

H x = - (cosh 2 77 - cos 2 gy 1 .j^sinh 77 * [cosh 77 + ( 1 - 2 /jl) cos - u(<f + ^s in £ • [( 

3 


(5) 


2^_ 


cosh 77 - cos £ 

1 - 2/jl) cosh 77 + cos - v(77)j ja 2 ^J. . 


L 



The functions u(£), v(^) are at our disposal. If they are taken as zero then h 0 is formally the 
same as the two-center Hamiltonian and consequently is separable. However, the more general 
H 0 given in Equation 5 is still separable and we have the functions u, v available to be chosen for 
our convenience , say to minimize certain perturbations, for example. 

Somewhat simpler expressions ultimately result if a contact transformation is made on H 
which effectively completes the square in the and polynomials inH 0 , the transformation 
being 


£ = 




> 


V 


A ! + 2 u(\j) , 

A 2 + 1 V { K 2 ) ■ 


V 


( 6 ) 


Then 


« = h 0 + n x . 


2 (cosh 2 k 2 - cos 2 /Vj) 1 Aj 2 + A 2 2 _ 4 u 2 (\J + ( 1 _ 2/-0 cos k x - 4 v 2 (\ 2 ) - cosh A 2 j 
- (cosh 2 k 2 ~ cos 2 1 ^sinh k 2 * [cosh k 2 + ( 1 ~ 2 fx) cos \ x ] - u(A t )| ^A 1 + ^ u (^1) 

+ |sin k x • [( 1 - 2 (i) cosh k 2 + cos A.J - v (A 2 )j- |a 2 + 2 v (^2) 


r < 7 ) 


j 


Thus the Hamilton Equation for the intermediate orbit is 


dt 


+ 2 (cosh 2 k 2 - cos 2 1 





+ ( 1 - 2 fi ) cos k x - 4 v 2 (\ 2 ) - cosh k 2 


0 * ( 8 ) 


By employing the separation of variables method we obtain as a complete integral 


~kt + 




2 (A) - ( 1 - 2 fx) cos A - ' 2 k cos 2 k - a dk + 


r*v5 


4 V 2 (A) + cosh k + 2 ^ cosh 2 A + a dA( 9 ) 


wherein A and a are canonical constants. 
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THE INTERMEDIATE ORBIT 


According to the Hamilton -Jacobi theory the intermediate orbit is provided from Equation 9 


by 


- Pi 

H 0 

— n P 

on 

-*2 

da 

A 1 


A 2 

" (9A.j 


> (10) 


where J3 19 p 2 are the remaining two canonical constants. 

In order to obtain a more specific class of orbits we shall specialize Equation 9 to some ex- 
tent by taking the functions u and v to be of the form 


4 u 2 (A.) 


4 v 2 (k) 


Uq + u 1 cos A + u 2 cos 2 k , 
v o + VjCoshA + V 2 cosh 2 k 


(ii) 


with u i , constants. This choice keeps the orbits within the general types that appear in the 
two-center problem. 

Introduce new constants i± i9 defined by 


Mo - u o “ a . 

Mi = Uj - (1 - 2fx) 


v n + cc , 


v \ ~ v l + 1 


M 2 


u 2 “ 2 A 


V 2 + 2 ^ 


> ( 12 ) 


From Equations 9 to 12 we then obtain the equations defining the intermediate orbit: 


* " /*! = 


cos 2 k dk 


Juoy 


±^0 + fj. t COS k + l± 2 cos 2 k 


cosh 2 k dk 




± + v \ cosh k + v 2 cosh 2 k 


(13a) 


5 



r a, 




d\ 


± |/p 0 + /Xj COS A. + p 2 COS 2 ^ 


d\ 


( 0 ) ± Vv 0 + v i cosh /V + v 2 cosh 2 


Aj - ± y/s 0 + p 1 cos +/x 2 cos 2 A.j 


A 2 - ± |/v 0 + v x cosh k 2 + v 2 cosh 2 X . 2 


(13b) 

(14a) 

(14b) 


We defer making the choice for the constants \. (0) in Equations 13a and 13b until we have ex- 
amined the character of the integrals more fully. 

The integrals are elliptic and cannot be expressed in compact form in terms of elementary 
functions apart from very exceptional values of the constants p i9 v.. The general character of 
the orbits depends most immediately on the zeros of the polynomials p 2 x 2 + x + p 0 and 
v 2 x 2 + Vjx + v 0 inasmuch as the vanishing of these polynomials provides the points of zero mo- 
menta by Equations 14a and 14b. The books by Charlier and Wintner (References 1 and 2) dis- 
cuss the dynamical significance of these roots in some detail for rather general problems. A 
proposed orbit for the Apollo capsule is frequently described rather loosely as being first an arc 
of a two-body ellipse extending from near the earth to the neutral point between the earth and the 
moon and then continuing on another two-body ellipse to near the moon, thereby describing a very 
slender figure of the form: j Now, orbits of this general character are included in Equation 13a 
and 13b: these are the cases i&x and 14/3 of Charlier (Reference 2, Vol. I, p. 126) which represent 
quasi-lemniscates winding about the two mass points. We shall assume 

m 2 > 0 
v 2 < o 



which is the case, for example, if u 2 , v 2 are zero or are quite small, and k < o. Let p l9 p 2 be the 
roots of fj . 2 x 2 + p 1 x +/x 0 = 0 and let r : , r 2 be the roots of v 2 x 2 + v x x + v 0 = 0. To give the orbits 
in question, have r 19 r 2 real with r x > 1 > r 2 and p 19 p 2 either complex or else real with p 1 >p 2 > 1. 
This means that the ranges of the variables \ x , k 2 are 0 < k x < 2v , o < cosh\ 2 < t x , since 
p 2 cos 2 k x + p x cos k l +p 0 = p 2 (cos k x -p x ) (cos k x - P 2 ) does not now vanish or go negative for real 
k 19 and v 2 cosh 2 k 2 + v x cosh X. 2 + v 0 remains positive or zero only in indicated range for k 2 . The 
explicit formulae for the roots in terms of the coefficients are 


~Pi + y M I ' 4 Po P 2 

2p 2 


Pi 4/x 0 p 2 

p 2 " 2/X 2 


~ v i - y^i 2 - 


0 "2 


2v, 


-*1 + - 4v o 

2v 0 


y (is) 
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It appears from the foregoing that a rather convenient choice for k^ and \ 2 <°> in Equations 
13a and 13b is \/ 0) = 0, \ 2 <°> = -r 1 where cosh r t - r x . The integrals in these equations thereupon 
vanish at the point (in barycentric coordinates) 


X 


~h L + 


1 I 

2 + 2 r i * 


y 


o , 


i.e., the point on the line connecting the two mass points, on the outside of the point with mass 
Thus p x is the instant of passage through this point. 


The solving of Equations 13a and 13b for k 19 k 2 in terms of t can be effected in terms of el- 
liptic functions but the expressions are prolix. Simpler formulae can be obtained through the in- 
troduction of parameters & 19 ^ 2 defined as 


dk 


1 V/ 2 0 + JjL x COS k + l± 2 cos 2 k 


rkr. 


d \ 


± 0 + v j cosh k + v 2 cosh 2 k 


(16) 


These are not independent, of course, since we have from Equations 13a and 13b, 


0-2 - ‘ 


The integrals in Equations 16 have relatively simple expressions in terms of inverse elliptic 
functions, and the equations solved for k x and k 2 give (Reference 4, p. 134): 


(<°2 + !) - 2 P 2 sr * 2 [2 " X ) ( p 2 + !) *1 

COS k x 1 

(p 2 + l) - 2 sn 2 >h 2 (pi - 1) (p 2 + *) A l- k l\ 

2r i ~ ( r i -1 ) sr > 2 [i V-2-2 ( r i " r 2 ) K ~^i)< ^ 2 ] 

cosh k 2 = r ^ ’ 

2 + (r t - 1) sn 2 || y-2v 2 (r t - r 2 ) (<►, -/3 2 ) , & 2 J 

r i h - r 2 ) + r 2 ( r i - i) 3 " 2 [ly _v 2 ( r i + x ) ( x - r 2 ) hi -^ 2 )* ^ 2 ] 
( i_r 2) + ( r i _i ) sr > 2 [W-M r » +i ) h _r 2 ) hi -^ 2 )' fe 2 ] 



2 hi zfA 

1 - !) h 2 + 




1 - l) ( r 2 + l) 

2 ( r 1 _ r 2 ) 


r 2 > — 1 - (17b) 



-l) (-1 t 2 ) 
+ 1 ) h- r 2) 


r 2 <- 1 . 

The function sn(u, k) is the Jacobi elliptic function, the sine amplitude of u with parameter k. The 
formula for cos Xj has all quantities pure real if we restrict ourselves to the case of p x , p 2 
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pure real, as we shall do from here on. Both formulae for cosh\ 2 are valid no matter 
whether r 2 >-lor r 2 <-l, but the form of the first is such that all quantities in it are real if 
r 2 > - 1, and similarly for the second if r 2 < - 1. Hereafter, though, we shall use only the first ex- 
pression for cosh k 2 . 

Equations 17a and 17b represent a pair of parametric equations with parameter for the 
intermediate orbit. If these expressions are introduced into Equations 14a and 14b formulae for 
the momenta in terms of result: 


A 


i 


A 


2 


2a 1 yp 2 ~ 1 dn (a : o x , k^ 
[p 2 + l) - 2sn 2 (a x k x ) ’ 


2a 2 yf 2 2 - 1 sn (a 2 » ^ 2 )^° ( a 2 ^2 * ^ 2 ) 

2 + ( r i - l)sn 2 (a 2 a 2 . k 2 ) 


a i = iy^2 C^i - !)(P2 + !) ■ 


a 2 - |V-2 v 2 (r, - r 2 ) , 


_ . 2 C°i Pi) 

1 " + ’ (18a) 


K = 


( r i ~ 0C r a + !) 
^ ( r 1 — r 2) 


b -/3 2 . (18b) 


By using the counterpart for elliptic functions of the trigonometric half-angle formulae 
(Reference 4, p. 146) we can express cos k x and cosh k 2 in a form devoid of squared elliptic 
functions: 


cos X.j 


cosh k 2 


-p 2 + 1 + 2p 2 cn (2a : h x } + ( p 2 + l) dn (2a x o x , k^ 

p 2 - 1 + 2cn (2a j + (p 2 + 1 j dn (2a x . k x j 


a i = " ! ) + ! ) ’ 


< 


0 _ / 2 ( p l l °2) 

1 = Wl~'){P2* 1) ’ 


(19a) 


r j + 1 + (r j - l) cn (2a 2 + 2r j dn (2a 2 <^ 2 5 ^2 ) 

r 2 + 1 - (r 2 - ljcn (2a 2 ^ 2 , & 2 ) + 2dn(2a 2 ^ 2 , & 2 ) 


r 



“ i )( r 2 + ~i)( r i T 2)- 1 


-^ 2 * 


Observe from Equations 17 or 19 that the degenerate case resulting from letting r x - 1 gives 
simply cosh k 2 - 1 which represents the straight line solution running from one of the mass points 
to the other. Also note that putting r x = 1 causes k 2 to vanish. And k x is small if p x is close to 
P 2 . Therefore, there are lemniscate-type orbits for which the elliptic function parameters k x and 
k 2 are small, which suggests expanding cos k X9 cosh k 29 A X9 A 2 in powers of k x and k 2 . The result 
is especially convenient since these expansions happen to involve elementary functions only. The 
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Fourier expansions of the elliptic functions (Reference 4, p. 147), 


V a n+1/2 

2K^sn(2Ku) - 477 / ^ "2^+1 s l n (2n + 1) * 

n - 0 


E a n 1/2 

l + ~ 2 nM cos (^ n + 1 ) 770 


2Kdn(2Ku) - 77 + 477 ^ ^ — 77 cos(2n77u) 

n = 1 


> ( 20 ) 


provide a means for obtaining representations valid for small k since 


«[■ 


k\ 6 


1 + 2I4J +15 tfi +150 lT>' + 


( 21 ) 


In Equation 20 K - K (k) is the "complete elliptic integral of the first kind": 


K<S> 5 = ““J- 


( 22 ) 


If q is replaced in Equation 20 through the use of Equation 21 there is obtained 
sn 2 Ku r sin77u + ^6 sin77u+ Y5 sin h 2 + o(& 4 ) , 
cn 2 Ku = cos 77u + cos 7 tu + Jg cos 37 tu^ h 2 + o(& 4 ) , 

dn 2 Ku - 1 + 4^ + 4" cos 27 tu^ k 2 + o(li 4 ) 


> (23) 


The formulae specifying the intermediate orbit valid for small k 19 k 2 are then found from Equa- 
tions 19, 18, and 20 to be 


COS 


cosh 


1 + p 2 cos Kj (p 2 2 - 1) 


77a t & 


1 11 


*1 


p 2 + cos £ 


3r l + 1 + ( r l - 1) 




4\p 2 + cos k< 
77a « 


V+o(v) , 


(24a) 


2 2 
COS £ 


( ri 2 - l) sin 2 


3 - (rj - 1) cos K 


2 c> 2 


( r i +3 ) " ( r i ~ !) cos 


k 2 + °{K) ’ (24t> ^ 


77a i 


A i ■ 2a i yPi 1- 


- — + 


(p 2 - l) + (p 2 - l) COS Kj 


p 2 +cos K 




4l p 2 + cos 


K, 


V r + 0 (V) 


(24c) 


A 2 - 4a 2 |/ ri 2 - 1 sin- 2 K 



~ 3r i “ 1 + ( 4r i +8 ) c° s “ K 2 "( r i - !) c 


k+3- (ri-l)cos K 


2 2 


r i+3-( r i~l) cos ' 


f, 2 


*>(V) • 


= (pi - 1) (p 2 + 1) • 


= V 2 (Pi -PjT^i " i)' 1 (^2 + i)’ 1 

K. = K(M , 

a 2 = if 2 v 2 ( r ;-r 2 ) , 


(24d) 


(24e) 


yr ( r i ■ oT 2 + 0 ( r i - r 2 ) _i 
*(* 2 ) , 


(24f) 


Avhere o - 2 = -/? 2 . 

Better (but more complicated) approximations than those of Equations 24a through 24f can be 
obtained by employing the process of "convergence improvement" on Equation 20. The success 
of the method depends essentially on the possibility of discovering a known Fourier series that is 
close to the given one. The difference series will then, in general, converge more rapidly than 
the original. We apply this method to the cosine amplitude functions. From Equation 20 we have 


Kk 

2tt cn(2Ku) - 


V -1 q n+1/2 

’ cos( 2n + 1 ) 77U 


L i 1 + q 

n = 0 


v r qn * i/2 

Z-. 1 + q 2n+1 

n = 0 L 

co 

L 


^n + 1/2 + qn+1/2 C OS(2n + l)77U , 


,n + l/2 


— * q3n+3/2 

" Z-i I + q 2n + 1 


cos( 2n + 1) 7 tu 


cos(2n + 1 ) 77u , 


( 1 - q) q 1/2 cos 77u 
(1 + q) 2 - 4qcos 2 77u 


-Ei 


3n+3/2 


c °s( 2 n + 1 ) 77 u 


( 1 - q) q 1/2 cos 77U 
( 1 + q) 2 - 4q cos 2 77U 


t!l v _ 

t( i Li 


.3n+3/2 


— COs( 2n + 1 ) 77U 


1 - q - q 2 + 2q 2 cos 2ttu - q 3 y~< q 3n+3/2 

q 7 cosTTU — 7TTTT/,' „ _i\ - / , , . 2rl +i cos(2n + 1 )ttu . 


(1 + q) (l - 2q cos 2rai + q 2 ) 


/ . 1 + q 3 "* 1 

n = 1 
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The Fourier series introduced above, 


L 


q n+l/2 


cos( 2n + 1 ) 7 tu 


was summed by employing simple operations on a known series (Reference 4, p. 213). A similar 
process works just as well for dnu so that we find 


K ft 


2tF cn(2Ku) 


1 /2 

q COS 77U 


- q - q 2 + 2q 2 cos 2mj - q 3 
( 1 + q) (l - 2q cos 2mi + q 2 j 


+ o{k 9 ) 


2K 

— dn(2Ku) 


(l “• 2q 2 ) + 2q (l - q 2 ) co s 2ttu + 
(l + 2q 2 ) - 2q (l + q 2 j cos 2ttu 




> (25) 


Consequently, these along with Equation 21 substituted into Equations 19a and 19b provide formu- 
lae for cos and cosh k 2 with errors of the eighth order in k 1 and k 2 . The resulting formulae 
are not written out here. 

The coordinates and momenta of the intermediate orbit are expressed in terms of the param- 
eter To locate position in orbit for a given instant it is necessary to know the relationship 
between and t. Under differentiation Equations 13a and 13b become 




0 


cos 2 A.j dk l 


cosh 2 k 2 d k 2 


t fio + Mi cos + (JL 2 cos 2 ± }f u o + cosh ^2 + y 2 cosh ^ 


d\. 


d k n 


± y^i 0 + cos k x + fi 2 cos 2 k 1 ± fao + v l cos b k 2 + V 2 cosh 2 k 2 


> 


J 


Also, from Equation 16, we find 


6k t 




+/lj COS \ j + fj. 2 cos 2 k y 


d & 2 - d-^ . 


Thus, 


d-t = i (cosh 2 \ 2 - cos 2 Xj) dd-j 


11 



or 


t - 0 , = 


cosh k r 


cos k. 


i r^i 

4 I (cosh 2 \ 2 - cos 2 du 

Jo 

2r l - ( r i _1 ) sn 2 a 2 (u-/3 2 ) 

2 + ( r l _1 ) sn2 a 2 ( u ~ @2) 

P 2 + 1 - 2p 2 sn 2 a : u 
p 2 + 1 - 2 sn 2 a x u 


> (26) 


There appears to be no very simple compact expression devoid of integrals for t in terms of . 
However, we can get the lead-off terms in an expansion in powers of k x and k 2 by using Equations 
24a through 24f and effecting the integration. The resulting formula is 


4 (*-/ 3 x ) 


P 2 (l ~P 2 f 


r 1 + 1 




1 {^fl ~ *) jp} ~ 1 


*7 - 1 + ^7 2 


k 2 

Tti 


arct an 


'P 2 ~ 1 'h *1 


P 2 + l 


tan ‘ 


2 KM 

/ r l + 1 2 ( r l + 1 ) 

7r 2 + 14r . + 27 

1 

to 

y 2 + ^2 c-i - 1 ) ■ 

4>^|^7TT ft 2 

^ 2 2 -i 

sin k x p 2 ~ 1 

sinft-j Oj ■ (2p 2 2 + 1 

Pi 

+ cos k x -&J + 4&-j 

(p 2 + cos k 

2(r7-l) 

sin^ 2 (<vi -/3 2 ) 

2 ( r i + i) ( r i- 

*2 • 

( r l + 3 ) _ ( r i _1 ) °° s M*i 

-Pi) Mi -1 )’ 

r 2 (n-i) 

/( r i+l) 2 ( r 1 + 1 ) 

1/ + P / i\ 1 

7r 2 + 14r. + 27 

— fcj 


arct an 


W 


tan 


Mv^) 


1 * 1 ; 




t[(fj - 1 ) cos^ 2 (^ - 0 2 ) - (r, + 3)] 2 




2 ^-l) 


2 ^2 ( r l *) 4^2 -j/rj + 1 

sinf>- 2 / 3 2 


arctan 


' r i +1 ^2 

o t an n 


\ ' ( r l +3 )“ ( r l _1 ) COS ' t 2^2 ' t 2( r l“ 1 )‘ 


2 (r, + l) ( r i “ !) sin ^2 ^2 •[ r i 2 + 2r i + 13 + (~ r i 2 " 6r i + 7 ) cos M2] 

4[( r 1 - 1 ) c os'& 2 /3 2 -( r i +3 )] 2 


M’"’ 


(27) 


where 


K N 


K(k 2 ) 
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The last two lines of Equation 27 represent a constant such that the entire expression vanishes 
with in accordance with Equation 26. Note that we have taken the constants of integration in 
such a way that k x = o and t = J3 t when = 0 . 

Equation 27 is not useful if k 2 is small because it is r t - 1 that is small, since the integrations 
have introduced r t - 1 in the denominators. This difficulty can be avoided in this case by express- 
ing the formula for cosh\ 2 in Equation 24b in terms of k 2 andn, 2 (instead of k 2 and ^ ) before sub- 
stituting in Equation 26 and integrating. 


PERTURBATION EQUATIONS 


The variation equations are formed on the Hamiltonian H l given in Equation 7: 

H x = - (cosh 2 \ 2 - cos 2 1 ^jsinh \ 2 • [cosh \ 2 +( 1 - 2p) cos k^ ^ u(\ 1 )| 

+ |sin • [( 1 - 2 fj.) cosh k 2 + cos \j] ~ v (^ 2 )j* j^ 2 + \ v (^ 2 

A little elliptic function manipulation on Equations 17a and 17b gives formulae for the trigono- 
■< ^ 
metric and hyperbolic sines which occur inH r : 



Similarly, 


sin Aj 


sinh k 


2 y P 2 ~ 1 sn a i cn a i 
P 2 ~ 1 + 2 cn 2 a t -a-j 

2 " ft t 2 ~ 1 cn a 2 & 2 

r, + 1 - ( r i - 1) cn 2 a 2 ^ 2 




u(A.j^ - ± y Uo + Uj cos X.j + u 2 cos 2 


+ + u i + u 2)^2 + l ) 2 - 2 (P 2 + 1) [2u 0 + u, (p 2 + 1) + 2 u 2 „] sn 2 a,, ^ + 4(u Q + u t P 2 +u 2 P 2 2 ) 


p 2 + 1 “ 2 sn 4 a t 


and 
l 


v ( x 2 ) = ± + Vj cosh k 2 + v 2 cosh 2 \ 2 


= ± 


V *( v 0 +V 1 r l + V 2 r l 2 ) + 2 ( r l - +v l( r l - 1)- 2v 2 r i] sn2 a 2 s + ( r l ~ 1 ) 2 ( V 0~ V 1 + V 2 ) ; 


*2 2 


2 + (r t - l) sn 2 a 2 o 2 


The constants u 4 , v £ are still at our disposal. For instance, the radicands inu(\,) and v(a. 2 ) are 
perfect squares if we take u 2 2 - 4u 0 u 2 = v t 2 - 4v 0 v 2 = 0 , and the formulae are somewhat simplified. 
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The equations of variation for the canonical constants are of the form 


dk 

(9H, 


d/3j 

<9H 1 

dt ~ 



“dt = 

- dk 



► 



da 

dHj 


d/3 2 

dH 1 

d-t 

’ 


d-t 

da 


(28) 


We may think of the coordinates and momenta as being expressed in terms of the parameter so 
that Hj in Equation 28 is to be regarded as known at this point as a function of ^ rather than of t. 
Consequently, we transform Equation 28 to get as the independent variable in place of t. Let 
us write k - a 19 a = a r Then, with the independent variable’s being t the variation equations are 


dcu 


A 2± 

at 


- IhT. »l[ A r (“i' l )' A r(“i' &V *)] 


> 


(29) 


The Hamiltonian for the problem we are treating is independent of any explicit value of K so we shall 
work out the transformation to of Equation 29 under this assumption. Let us represent the re- 
lation between ^ and t by 


t - <p(a i , 

Oj - 0(a it /? it t.) • 


(30) 


In our particular case this transformation is, of course, that of Equation 26. The total differential 
of t is, according to Equation 30, 




d<p 

da 


da 


d<p 

M r 


d/? 


d<p 

+ d ^l 


in which the summation convention is observed, i.e., a term having an index repeated is to be re- 
garded as a sum taken over that index. Thus, the variational equations take on the form 


da . 

da 

r 

(9Hj 

da-j 

^ i r d^j 

= 

d^ 


<9H 

d-a 1 

^ i r da- j 

da 


dcfc da r d0 ^ t d& 

da d-a-. + d(3 d-a. + d&* 

r 1 ' r 1 1 


d<t> da r d& 9<t> 

da r dOj + d/3 T d-&j + 
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in which S. r is the Kronecker delta: S. r = 0, i / r, and S ir = l, i = r. In a standardized form 
the equations are 


<9H 


I d4> 


da 


r d/3, 

dfi i da r - s ir) + <9/3. 5/3 r dOj 


c9H 


1 <90 


da r 


<9 a i <?a f do-j 


j90 

+ (9a. <9/3 r " 8 ir) d^j 


<90 

<9/3. <9^ 


Mi <90 
<9a. (9^ 1 


(31) 


where 


a x - ^ 


a 9 cl , 


wherein the summation convention applies. It remains to be seen how the various quantities in 
Equation 31 are to be computed. Now 


H x - H 


l [M a i- t )’ A r ( a i' *)] " Hi [\ r ( ai , /3., 0(a., /3. , A,)) , A r (a . , /3. , </> (a ; , /3. , 

s Hj (a r , /3 r , Oj) = Hj (a r> /3 r , >/((a r , / fl r , 


and it is this function n x (a r , j3 r , ) that the Hj in Equation 7 actually represents. The <9H l /dj3 i and 
dH x /da. that occur in the coefficients in Equation 31 are thus found from 


(32) 


<9H X 

<9H 1 

<9Hj 

<90 


i ^ 

ii 



<9H 1 

<9Hj 

<9H 1 

<90 


<9a . 

i 

+ <9^j 

<9a. 


with t then replaced by 0(a r , /3 r , ) . As a matter of fact, <90/<9/3. and <90/<9a. will be computed in 

practice by implicit differentiation of K = 0(a., /3. , where s 0, which automatically replaces 
* by (a function of) * t . The derivatives d<p/da r , <90/<9/3 r , <90/<9^ are simply the straight partial de- 
rivatives of the function 0(a. , /3., . 

We shall not carry out the details of determining explicitly the coefficients in Equation 31 for 
our particular case but shall content ourselves with working out necessary auxiliary formulae. 

We find that 



-> 

1 


dv 2 

1 

dti ~ 

2 ’ 

> 

dh ~ 

2 

£Mo 

- 1 , 
y 


d Z o 


<9a 


da 

1 




> 


(33) 
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with all remaining derivatives of p, i 


v i with respect to a, k, and J3. zero. Furthermore, 


dp. 

Pi 


dp i 

1 

1 

r 

i dr, 1 , 

dk 

4 ’ 


da 

2a j 2 k 2 

dA 

4a 

2 <?a 9 2 

2 ^ a 2 

dp 2 

pi 


dp 2 

1 

ar 2 


r 2 2 5r 2 1 

dk 

4a j 2 k 2 


da 

' 2 3l 2 V 

a-fi 


4a 2 2 ’ da " 2a 2 2 





*Pi 

dr. 







M, = 

dH, - 

0 , 




da 1 

pi - 

- p 2 + 2pj p 2 


<9 a 2 

F 1 r 2 



dk 


16aj 


<9fi 

4a 2 ( r 1 — r 2 ) ' 



da i 

2 - 

Pi. + ^2 


da 2 

1 



da 

8a l( 

>1 "Pa) 


<9a 

2 a 2 (r , ~ r 2 ) 





<9 a . 

i 









= 0 , 




dk^ 


{-pi -p 

2 2 + 2 pi pi) Aj 



r, 2 + r 2 2 - 2 Tj 2 r 2 2 


dk 


168^ 



dk 

163 2 2 (r, -r 2 ) 2 ^ 2 


dk l 


[pi + pi 

2)^1 


dk 2 

r 2 2 + r 2 2 - 2 


da 


8a i 2 (pi 

-^) 2 


da 

8a 2 2 ( r l " r 2) 2 k 2 


dk. 

= 0 • 





d cosh k 2 

~~ M. 


[ 2 + ( r » " 1) sn2a a( u -^rfe- J a 2 (u-/i 2 ) " 


<9sna 2 (u-/3 ) dlt 9 l 

^2) cn a 2 (u-^)dn a 2 (u -^) - 4 (r , 2 - l) sn a 2 (u - ^j) gk | 


d cosh A. 


da 2 - [ 2 + ( r i- i) sn 2 a 2 (u-/3 2 )] 


-2 2 


2 (rj 2 - l) 

i( u "/ e 2 ) " a 2 (rj _r 2 ) ( u “^2) sna 2 ( u 


d COS k ^ 
<9fi 


. d sn a (u “ j3 0 ) 

/3 2 ) cna 2 (u-/3 2 ) dna 2 (u-/3 2 ) - 4 ( r j 2 -l) sna 2 (u-/3 2 ) gh , 


9 “|— 2 | P 2 r\ 

2 sn a x uj < 2 ^ 2 sn2 a i u cn 2 aj u 


[>2 + 1) 

W - 1) (p, -P 2 + 2 Pl p 2 ) 


"1 U 


4a . 


<9 sn 3j u dh j j 

! i ucnaj udna, u-4(p 2 - l) sn a, u ’ "jF , 


( 41 ) 


<9 co s A. 
<9a 


j^(p 2 + 1 ) ~2sn 2 a 1 uj 


2 2 

2 sn^ a, u| < 2 sn 2 a 1 ucn 2 a 1 u 

a i N 


t°2 0 (^i P2 + P 2) d sn a x u (9&J 

2a i(/°i "P 2 ) U Snai UCna i udna i u '4(p 2 2 - l) snaj u ^ 


CONCLUDING REMARKS 

The formulae of the preceding section are lengthy and involved, and it appears difficult to get 
pertinent information out of Equation 31. For a first order theory it may well be better to deal 
with the variation equations in t of Equation 28, rather than those in of Equation 31. The useful- 
ness of the partitioning introduced in this paper will, of course, depend largely on how much trouble 
it is to handle the perturbation equations and on how large the perturbations turn out to be. As far 
as application to the Apollo project is concerned, work remains to be done to ascertain the prob- 
able range of values for r i? p. for these orbits, so that it can be determined whether the type of 
approximation used in representing the intermediate orbit (i.e., low order in and h 2 ) will be 
reasonable. 


A certain amount of latitude is present due to the arbitrary character of the functions u(\j) 
and v(\ 2 ) introduced at the partitioning of the original Hamiltonian. This paper considers a spe- 
cific choice for the form of these functions, although six arbitrary constants are left available. 

An investigation could well be made concerning the feasibility and usefulness of some other choice 
for these functions. 


17 



REFERENCES 


1. Wintner, A., "The Analytical Foundations of Celestial Mechanics," Princeton, N. J.: Princeton 
University Press, 1947. 

2. Charlier, C. V., "Die Mechanik des Himmels," Leipzig: Veit and Co., 1902. 

3. Samter, H., "Das Zweizentren- Problem in der Storungstheorie," Astron. Nachr . 217(5193- 
5194):129-152, April 10, 1922. 

4. Magnus, von W., and Oberhettinger, F., "Formeln und Satze fur die Speziellen Funktionen der 
Mathematischen Physik," Berlin: Springer-Verlag, 1948. 


18 





Appendix A 


List of Symbols 

Two finite masses in the restricted three-body problem 

x, y Barycentric coordinates for a system rotating with masses such that particle coordinates 
for fj. and 1 - jj. are (1 - 0) and {-jj., 0) 

L, u Lagrangian and force functions for the restricted three-body problem 

r j , r 2 Bipolar coordinates representing distances of (x, y) from the points (- m, 0) and (1 - fi, 0), 
respectively 

Hi r) Elliptic coordinates expressed in terms of r l9 r 2 
f £ Momenta conjugate to r. 

a Canonical constant in Equation 9 

, j3 2 Canonical constants in Equation 10 
Mi, tq Constants defined by Equation 12 
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